This document summarizes updated results from the analyses described here. Notably, the analysis now includes Pink (n=70) and Chum (n=44) salmon stocks from across the eastern North Pacific, as well as updated data for Sockeye stocks (n = 52).
The overarching goal of our analyses is to characterize relationships between ocean conditions (SST and abundance of potential competitors) and salmon productivity (R/S) over space and time. Since our last update, we have made minor changes to our approach, as described in the sections below.
We have compiled spawner-recruitment time series for all five species of Pacific salmon, which are hosted in a github repository, although only Sockeye, Pink, and Chum populations are analyzed here. We analyze a total of 166 populations originating from Washington to the Bering Sea:
From our original set of 58 Sockeye populations, we obtained more data for the Southeast Alaska and West Coast (BC) regions, but subsequently chose to remove time series that did not span at least 1984-2015 brood years without significant gaps. This is intended to make sure each time series overlaps with all three time periods in our era models, and excludes short series that random walk models may be sensitive to. There are now 52 Sockeye stocks in the analysis - many of the original stocks, and some new ones.
We obtained spawner-recruitment time series for Pink and Chum populations, and currently use all time series with good data quality, regardless of whether they meet the criteria applied to Sockeye data. These data are not representative of all regions and time periods, and notably, we have no data after 2009 for any Alaskan stocks. We are working to update these data through to the present.
Compilation of coho and Chinook time series is ongoing.
The map below shows the ocean entry locations of all stocks in the analysis presented. The legend to the right of the plot allows you to toggle between species (double-click a legend item to isolate that species). Hovering over a point will show the stock name and ocean region.
Figure 1. Ocean entry locations of Sockeye (n=52), Pink (n=70), and Chum (n=44) salmon stocks included in analyses (total N=166).
Productivity (log[R/S]) time series of all stocks are shown below to illustrate the length and relative number of time series among regions and species. Vertical dashed lines represent proposed regime shifts in the ‘era’ model (described in sec. 2.2) to illustrate the coverage of data over periods of interest. Where coloured lines break to light gray, data are missing.
Figure 2. Productivity (log[R/S]) time series of Sockeye stocks (n=52) with vertical dashed lines indicating proposed ocean regime shifts.
Figure 3. Productivity (log[R/S]) time series of Pink stocks (n=70) with vertical dashed lines indicating proposed ocean regime shifts. Note that Even- and Odd-run pinks are considered separately due to their independent and sometimes vastly different productivity regimes.
Figure 4. Productivity (log[R/S]) time series of Chum stocks (n=44) with vertical dashed lines indicating proposed ocean regime shifts. Note that there are no time series in the Southeast Alaska (SEAK) region; one stock (Middle Skeena) that would have been assigned to SEAK was grouped with West Coast stocks instead.
Our updated analyses focus on three classes of generalized spawner-recruitment models:
Stationary models (e.g. Connors et al. 2020), which estimate time-invariant relationships. These have not changed since last update other than the addition of data for the two new species.
‘Era’ models (e.g., Malick 2020), which allow relationships to vary among pre-defined periods that represent hypothesized shifts in North Pacific Ocean processes and relationships.
We no longer model 1976/77 as a regime shift as per feedback from our last update, which highlighted that this shift from a cold to warm PDO phase simply reflected a temperature change (which is already accounted for via our direct consideration of ocean temperature) and not a fundamental change in oceanographic processes like the 1988/89 regime shift.
Our updated analyses consider a potential regime shift at the onset of the ~2013-2015 marine heatwave (‘the blob’). Currently, we consider brood years >= 2011 to have interacted with the marine heatwave regardless of the species and stock. In future analyses, the timing can be altered to better reflect the diverse life histories of populations in the analysis.
Random walk models (e.g., Malick 2020), which allow relationships to evolve gradually through time.
Details on each model class and our Bayesian model fitting procedure can be found here.
Across the three species, we continue to see evidence for spatially and temporally variable responses to warming and competition at sea.
In Sockeye, we continue to see strong spatial stratification of relationships with both SST and competition, with the southernmost stocks initially having the most negative responses to temperatures (and northernmost stocks having positive relationships). However, this structure disappears in the ‘marine heatwave’ period, in which Alaskan (SEAK and GoA) stocks have more negative SST relationships than more southern (WC) stocks, and all only WC stocks show negative responses to competition. However, due to the short duration of the marine heatwave period, we should interpret these findings cautiously.
In Pink salmon, there are clear spatial and temporal patterns in SST effects, but less pronounced patterns in competition effects, which are either uncertain (in early ‘era’) or near zero (in middle era). This is consistent with the hypothesis that Pink salmon are the least affected by competitors at sea (Ruggerone & Nielsen, 2005).
..chum insights to come with final model runs…
Figure 5. Posterior probability distributions of the predicted effect of SST (top), competitor abundance (middle), and the combined effect (bottom) on Sockeye productivity (R/S). Faint lines show stock-specific effects while bold lines show regional effects from hierarchical model. X-axis values represent the percent change in productivity per standard deviation unit increase in the covariate.
Figure 6. Posterior probability distributions of the predicted effect of SST (top), competitor abundance (middle), and the combined effect (bottom) on Pink salmon productivity (R/S). Faint lines show stock-specific effects while bold lines show regional effects from hierarchical model. X-axis values represent the percent change in productivity per standard deviation unit increase in the covariate.
Figure 7. Posterior probability distributions of the predicted effect of SST (top), competitor abundance (middle), and the combined effect (bottom) on Chum salmon productivity (R/S). Faint lines show stock-specific effects while bold lines show regional effects from hierarchical model. X-axis values represent the percent change in productivity per standard deviation unit increase in the covariate.
Figure 8. Posterior probability distributions of the predicted effect of SST and competitors on Sockeye productivity over three pre-defined time periods/eras (earliest in top panel). Regional mean effects are shown by bold lines and individual stocks’ distributions by light lines.
Figure 9. Posterior probability distributions of the predicted effect of SST and competitors on Pink salmon productivity over three pre-defined time periods/eras (earliest in top panel). Regional mean effects are shown by bold lines and individual stocks’ distributions by light lines. Note that data from Alaska for the most recent time period are needed (Figure 3).
Figure 10. Posterior probability distributions of the predicted effect of SST and competitors on Chum salmon productivity over three pre-defined time periods/eras (earliest in top panel). Regional mean effects are shown by bold lines and individual stocks’ distributions by light lines. Note that data from Alaska for the most recent time period are needed, and there are no time series in the Southeast Alaska region.
Figure 11. Time-varying posterior mean estimates of SST and Competitor covariate effects on Sockeye productivity, modelled as a random walk. Individual stock estimates are in faint lines, while regional means are represented by bolded lines. Region-wide means are post-hoc calculations, rather than resulting from hierarchical model structures as in the stationary and era models.
Figure 12. Time-varying posterior mean estimates of SST and Competitor covariate effects on Pink salmon productivity, modelled as a random walk. Individual stock estimates are in faint lines, while regional means and 80% CI are represented by bold lines and shaded areas. Solid bold lines represent odd-year Pink stocks, while dashed lines are even-year stocks. Region-wide means are post-hoc calculations, rather than resulting from hierarchical model structures as in the stationary and era models.
Figure 13. Time-varying posterior mean estimates of SST and Competitor covariate effects on Chum salmon productivity, modelled as a random walk. Individual stock estimates are in faint lines, while regional means are represented by bolded lines. Region-wide means are post-hoc calculations, rather than resulting from hierarchical model structures as in the stationary and era models.
Priorities for future analysis are:
Update data to present (as recent as possible) for all species and regions, particularly Alaskan pink and chum
Sensitivity analyses of SST and competitor index covariates
Revise hierachical grouping of Pink salmon stocks | The random walk models for Pink salmon show the trends for even-year and odd-year stocks separately for ease of visualization, but the difference in the effects between lineages in some cases begs the question of whether even and odd stocks within a region should be grouped together in hierarchical models or separately. We will investigate the possibility of separating even and odd lineages, which means treating them as if they were separate species.
Random Walk model improvements | Work on versions of random walk models with (1) estimated region-wide mean effects, and (2) reduced sensitivity to time series of different lengths, is ongoing.